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Abstract 

We present a simplified version of the DFS-based Left-Right planarity 
testing and embedding algorithm implemented in Pigale ^12, which has 
been considered as the fastest implemented one 0. We give here a full jus- 
tification of the algorithm, based on a topological properties of Tremaux 
trees. 

1 Introduction 

It is well known since the publications 0J|3] in 1973-74 by J. Hopcroft and R.E. 
Tarjan that the time complexity of the problem of graph planarity testing is 
linear in the number of edges. However, authors of text books and teachers 
of graph theory know how hard it is to describe and completely justify such 
algorithms, each one being more tricky than the other. To display a planar 
embedding is usually a second hard step UJ. 

A satisfactory graph planarity algorithm should be efficient, clearly justified, 
and easy to understand. Only a strong mathematical insight of the subject could 
probably meet these three criteria. 

R. Tarjan initiated the use of the Depth-First Search procedure (DFS) to 
attack planarity testing, both DFS and planarity testing on biconnected graphs 
being handled recursively. Several authors explained how to extract an embed- 
ding during the recursion of such an algorithm. The authors produced in the 
eighties the so-called Left-Right algorithm, a non-recursive version avoiding 
the 2-connectivity assumption. The Left-Right algorithm appeared to be 
extremely efficient for testing planarity and embedding planar graphs and has 
been recognized as the fastest among the implemented ones by the comparative 
tests performed by graph drawing specialists But even if a tentative math- 
ematical explanation appeared in our papers on Tremaux characterization of 
planarity 00 El, it was not enough to fully meet the second and third criteria. 



The code appears in the GPL- licenced software PIGALE ^ 12 ■ A new, sim- 
plified, and faster version has been implemented at the occasion of this paper, 
which is the one we will discuss here. 

In Section [21 Tremaux trees are studied as a mathematical object per se. A 
rooted tree that spans a graph is a Tremaux tree if each cotree edge is incident to 
two comparable vertices (with respect to the tree order). It is why cotree edges 
are then called back-edges. For a given Tremaux tree, the following structural 
concepts are defined: the low of a vertex or of an edge, the low set of a vertex 
or of an edge and the fringe of an edge. 

Section 13 is devoted to the study of planarity. Planarity has been related to 
Tremaux trees in in a characterization based on the existence of a bicoloration 
constrained by three special simple configurations. These configurations are 
here unified into a single one. In this setting, it is proved that some additional 
constraints may be imposed, which do not change the existence of a bipartition 
but lead a simple planar embedding. 

The planarity testing and embedding algorithm is then described in Sec- 
tion 0] It is shown that linear time-complexity is reached by implicitly building 
a spanning arborescence in the graph of constraints. It shall be noticed that 
the data structures used by the algorithm are almost trivial ones, which may 
explain its computational efficiency and the ease to follow the algorithm step 
by step on an example. 

In the sequel, by a "graph" we mean a "connected, loopless multigraph", 
unless we state it otherwise. 

2 Tremaux Trees 

Depth-first search (DFS) is a fundamental graph searching technique known 
since the 19th century (see for instance Luca's report on Tremaux's work |10j) 
and popularized by Hopcroft and Tarjan m the seventies. The structure 

of DFS enables efficient algorithms for many other graph problems [U^. Per- 
forming a DFS on a graph defines a spanning tree with specific properties (also 
known as a Tremaux tree) and an embedding of it as a rooted planar tree, the 
edges going out of a vertex being circularly ordered according to the discover 
order of the DFS). 

A rooted spanning tree T of a graph G defines a partition of the edge set 
of G into two classes, the set of tree edges E(T) and the set of cotree-edges 
E{G) \ E(T). It also defines a partial order ^ on V(G): x -< y if the tree path 
linking y to the root of T includes x. The rooted tree T is a Tremaux tree if 
every cotree edge is incident to two comparable vertices (with respect to -<). A 
Tremaux tree T defines an orientation of the edges of the graph: an edge {x, y} 
(with x -< y) is oriented from x to y (upwards) if it is a tree edge and from y to 
x (downwards) if it is a cotree edge. Cotree edges of a Tremaux tree are called 
back edges. We will denote by lo + (v) the set of the edges incident to a vertex v 
and going out of v. When T is a Tremaux tree, the partial order < is extended 
to V(G) U E(G) (or to G for short) as follows: for any edge e = (x, y) oriented 
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from x to y, put x ~< e and if x -< y (that is: if e is a tree-edge) also put e <y. 

(6,1) (6,3) 




Figure 1: The partial order -< defined by a Tremaux tree of K 3j3 (cover relation 
correspond to bottom-up arrows as usual). 

Notice that in the partial order -<, all elements a and (3 of G have an unique 
greatest lower bound (meet) a A (3. Moreover, the maximal chains of -< all have 
the same structure: they begin with the root vertex of T, alternate between ver- 
tices and edges, and include at most one back-edge (which is then the maximum 
of the chain). 

When a -< (3, the unique chain (of -<) with minimum a and maximum (3 
which is maximal (with respect to set-inclusion) is denoted by [a ; (3] . 

It will be helpful to introduce a notation for the minimal element of the 
interval ]a ; 0\ : 

Definition 2.1 For x <e, where x E V(G) and e e E(G) we define 

stem(a;, e) = min ]x ; e] (1) 

This means that f is the first edge in the unique chain of ~< with minimum x 
and maximum e. 

Notice that in this definition, as in the remaining of the paper, intervals and 
"min" will always be related to the partial order -< defined by the considered 
Tremaux tree. 

Definition 2.2 The function low : G — > V(G) is defined by 

low(a) = min ({a} U{u£ V(G) : B(u, v) e E(G) \ E(T), (u, v)ha}). 
Notice that this function is well defined and that low(a) ^ a. 
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In the literature, this function is usually only defined on V(G). Our extension 
is such that for any edge e = (x, y): 

, . J min(x, \ow(y)) if e is a tree edge, 
\ y if e is a back-edge. 

Definition 2.3 The fringe Fringe(e) of an edge e = (x,y) is defined by: 

Fringe(e) = {/ G E{G) \ E(T) : f h e and low(/) ~< x } 

Definition 2.4 The low set L(a) of a e V(G) U E(G) is defined by: 

L(a) = {/ e E(G) \ E(T) : f h a and low(/) = low(a) } 

Notice that L(e) = {e} if e is a back-edge. 

Definition 2.5 Given a tree edge e e E(T), we 
will call e = (x, y): 

• a block edge 

if x < low(y), that is if low (e) = x (this 
means e is either an isthmus or the mini- 
mum edge of a block of the graph); 

• a thin edge 

i/low(y) -<, x and there exists no back-edge 
(u, v) with u >z y and low(y) -<u < x; 

• a thick edge 

if\ow(y) ~< x and there exists a back-edge 
(u, v) with y and low(y) ~< u -< x. 

Definition 2.6 A TT-precedence order -<* is a partial order on E(G) such 
that, for any v G V(G) and any e, / G uj + (v): 

• if low(e) -< low(/) then e ~<* f, 

• i/low(e) = low(/), f is a thick tree edge but e is not, then e ~<* f. 
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Fig. 2. Block, thin and thick edges. 

Block edges are dotted, thin 
edges are light and thick edges 
are fat. 
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3 Tremaux Trees and Planar ity 



Planarity has been related to Tremaux trees by de Fraysseix and Rosenstiehl 
in a series of articles |5] ■ One of these characterizations is based on a the 
existence of a special bipartition of the low angles of the back-edges into left 
ones and right ones. The constraints that the bipartition has to fulfill is encoded 
into two relations, namely the T-alike and T-opposite relations. 



O vertex 
| tree edge 

back-edge 

non trivial tree chain (at least one edge) 

possibly empty tree-chain 
< (starting with an edge, ending with a vertex) 
> 




subtree possibly reduced to a vertex 
tree path possibly reduced to a vertex 



Figure 3: Symbols used in the figures 

We don't give here the formal definition of T-alike and T-opposite relations 
in terms of -<, but simply recall the characterization given in [xj: 

Theorem 1 Let G be a graph with Tremaux tree T . Then G is planar if and 
only if there exists a partition of the back-edges of G into two classes so that any 
two edges belong to a same class if they are T-alike and any two edges belong to 
different classes if they are T-opposite. 

Instead of working with this characterization, we introduce an equivalent 
characterization based on a single configuration. 

Definition 3.1 Let v be a vertex and let e±,e2 G lu + (v). 
The interlace set Interlaced(ei, e2) is defined by: 

Interlaced(ei, e^) — {/ G Fringe(ei) : / y low(e2) } 

Definition 3.2 Given a graph G and a Tremaux tree T of G, a coloring A : 
E(G) \ E{T) — > { — 1, 1} is an F-coloring if, for every vertex v and any edges 
ei,e2 G uj + {v), Interlaced(ei, e2) and Interlaced(e2, ei) are monochromatic and 
colored differently. 
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case (i) case (ii) case (iii) 

a and f3 are T-alikc a and (3 are T-oppositc a and (3 are T-opposite 



Figure 4: Definition of T-alike and T-opposite relations 

It is easily checked that a coloring is an T-coloring if and only if two T-alike 
back-edges are colored the same and two T-opposite back-edges are colored 
differently. 

In a planar drawing, an T-coloring is defined by the partition of the back- 
edges on two sets, the edges / having their low incidence on the left (resp. the 
right) of the tree edge stem(low(/), /). The following lemma is straightforward 
and does not deserve a proof: 

Lemma 1 Let G be a planar graph with Tremaux tree T. Then G has an F- 
coloring. □ 

Definition 3.3 An F-coloring A : E(G) \ E(T) — > {— 1, 1} is strong if, for any 
v £ V(G) the low set L(v) is monochromatic. 

Lemma 2 If G has an F-coloring then it has a strong F-coloring. 

The addition of the constraints that the sets L(v) are monochromatic may not 
lead to a contradiction, as the only constraints involving f\ £ L(t>) would also 
involve any fi £ L(w) and would require that f\ and fi actually have the same 
color. □ 
^From the "low angles bicoloration" A we define a "high angles bicoloration" 
A on the whole edge set of the graph (on both tree edges and back-edges). 
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Definition 3.4 Let A : E(G) \ E(T) -> {-1,1} be a strong F-coloring. We 
define the coloring A : E{G) — ► {— 1, 1} by: 

c-, s J A(e), if e is a back-edge 

\ ^(/)) if e is a tree edge and f G Fringe(e) with maximal low(/) 

Lemma 3 Let G be a graph, let T be a Tremaux tree of G, let X be a strong 
F-coloring and let A be the associated mapping. We define the circular order of 
the edges at a vertex v as follows: 

Let ei 62 y* ■ ■ ■ e p be the edges in lo + (v) with A(e^) = —1 and let 
e p+i ^* e p+2 ~<* ■ ■ ■ e q be the edges in lo + (v) with A(ej) = 1. 

In the circular order around v one finds the incoming tree edge (if v =/= r) 
and then Li, e\, R\, L 2 , e 2 , R 2 , ■ ■ ■ , L q , e q , R q where Li (resp. Ri) is the set of 
incoming cotree edges f = (x,v) such that X(f) = —1 (resp. A(/) = 1 ) and 
the tree-path linking r to x includes e^. For ei,ej G L^ (resp. for e^e^ G Rk), 
one finds before (resp. after) ej in the circular order if stem(e,A, ej, ej) ~<* 
stem(ei A ej, ej). 

Then these circular orders define a planar embedding of G. 

In a drawing where the tree edges cross no other edges, only two kind of 
crossings could occur: 




Let x = e A f,y = low(/),z = low(e),a = stem(x, e) and (3 = stem(x,/). 
Without loss of generality we assume A(/) = — 1. 

• The first case corresponds to two back-edges e, / with A(e) = A(/) and 
low(e) < low(/). 

If low(a) = low(/3) = low(e) = low(/) = u, then the contradiction 
arises from the definition of the circular order at u. Otherwise, as / G 
Interlaced(/3, a), / is colored the same way as the back-edge in Fringe(/3) 
with maximal low-value. Hence A(/3) = A(/) = —1. According to the def- 
inition of the circular order at x, A(a)) = — 1 and a >-* (3. According to 
the thin-thick precedence of <* it follows that the edge e' G Fringe(a) 
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with maximal low-value is such that low(e') >~ low(/3), that is: e! £ 
Inter laced (a, /?). Thus A(a) = A(e') ^ A(/), a contradiction. 

• The second case corresponds to two back-edges e, / with A(e) ^ A(/). By 
symmetry we may assume low(a) < low(/3) and A(e) = 1. 

Then low(/3) >- low(x) as A is a strong F-coloring. Thus / is colored 
the same way as the back-edge in Fringe(/3) with maximal low-value and 
A(/3) = A(/) = —1. According to the definition of the circular order at x, 
A(a)) = —1 and a >~* (3. It follows that e e L(x), a contradiction. 

□ 

Theorem 2 Let G be a graph, let T be a Tremaux tree of G. The following 
conditions are equivalent: 

(i) G is planar, 

(ii) G admits an F-coloring, 

(iii) G admits a strong F-coloring. 

Moreover, if G is planar, any strong F-coloring A defines a planar embedding 
of G in which a back-edge e has its lowest incidence to the left of the tree if 
A(e) = —1 and to the right of the tree if A(e) = 1 

This is a direct consequence of the previous lemmas. □ 

4 The Planarity Testing Algorithm 
4.1 Outline 

Let G be a graph of size m. The three steps are performed in 0(m)-time. The 
first step is composed of a preliminary DFS on G and the computation of the 
low function and the status of the edges (block/thin/thick). The second step is 
the computation of a TT-precedence order, which may be efficiently performed 
using a bucket sort. We now examine the last step of the algorithm, which tests 
the planarity of the graph. 

We shall consider some data structure CS responsible for maintaining a set 
of bicoloration constraints on a set of back-edges. We assign to each edge of the 
graph e such a data structure CS(e). These structures are initialized as follows: 
CS(e) is empty if e is a tree edge and includes e (with no bicoloration constraints) 
if e is a back-edge. We say that all the back-edges have been processed and that 
the tree edges are still unprocessed. 

• While there exists a vertex v, different from the root, such that all the 
edges in u + (v) have been processed. Let e = (u,v) be the tree edge 
entering v. Let e\ -<* e 2 . . . -<* be the edges in uj + (v) (k > 1). We 
do the following: 
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- Initialize CS(e) with CS(ei). 

- For i : 2 -> k, Merge CS(e,) into CS(e), that is: add to CS(e) the 
edges in CS(ej) and add the F-coloring constraints corresponding 
to the pairs of edges ej , ej with j < i (notice that all the concerned 
back-edges belong to CS(e)). If some constraint may not be satisfied, 
the graph is declared non-planar. 

— Remove from all of the CS(e) every back-edge with lower incidence 
u. 

— We declare that edge e has been processed. 

• As all the edges have been processed, we declare that the graph is planar. 

4.2 Data-structure, Complexity and Embedding Compu- 
tation 

First notice that the processing order of the tree edges is simple to compute, 
using either a "topological sort technique" (by maintaining the unprocessed 
outdegree of vertices) or by following the backtrack order of a DFS steered by 
the TT-order. Also, we should notice that the structure CS(/) is only used 
(after its computation) when computing CS(e) where / is the predecessor of / 
in -<. It follows that we may "destroy" CS(/) when computing CS(e) without 
any risks. 

The data structure for CS(e) is a stack of double top-to-bottom linked stacks. 
This means that each element of CS(e) is a pair (S^S 1 ) of stacks. Each such 
pair correspond to a complete bipartite constraint graph: all the edges in S* (i E 
{0, 1}) have to be colored the same and they have to be colored differently from 
the back-edges in S J (for j ^ i). These will be all the constraints encoded by 
the CS(e) (hence no constraint exists between back-edges belonging to different 
stack pairs) . 

Some additional constraints will allow to get amortized constant time oper- 
ations: Let (S 1 ?, S{), (5§, S%), . . . , (Sg, Si) be the pairs stacked in some CS(e). 
Then: 

• The top back-edges of Sf of Sj have both a lower incidence which is strictly 
smaller than the lower incidences of the bottoms of Sf +1 and Sh^. 

• the lower incidences of the back-edges belonging to some stack Sf are in 
non-decreasing order. 

Notice that, for any tree-edge e = (u, v), CS(e) will exactly include the back- 
edges / = (x,y) such that x >- v and u >~ y. Hence the set of edges in CS(e) is 
Fringe(e). 

Now consider the operations performed by the algorithm on the structures 
CS(e): 

• Merge of CS(ej) into CS(e): the structure CS(e) then contains the back- 
edges in X{e{) — Uj<i Fringe (e^) and their bicoloration constraints. By 
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induction, the graph of these constraints is a disjoint union of complete 
bipartite graphs encoded by the pairs of stacks. The F-coloring conditions 
now express as: 

— all the back-edges in CS(e) (i.e. in X(ei)) which have a lower inci- 
dence strictly greater than low(ej) should have the same color, 

— all the back-edges in CS(ei) which have a lower incidence strictly 
greater than low(u) should have the same color, 

— the two above sets of edges should be colored differently. 
These constraints are added as follows: 

— let a be the smallest lower incidence of the back-edges in CS(e) (ac- 
cording to monotonicity we just have to look at the bottom back- 
edges of the stacks of the bottom pair). Then a = low(u). 

— let b be the smallest lower incidence of the back-edges in CS(ej) (as 
above, this is computed in constant time). Then b = low(ej). 

— if CS(e) = ((Sf, Si), (Si S^,...,(S° k ,Sl)), let j be the biggest in- 
teger < k such that none of the back-edges in the top of S® and Sj 
have a lower incidence greater than b. Then it should be checked 
that for j + 1 < j' < k, one of Sj, and Sj, is empty (for other- 
wise, we have found a contradiction in the constraints proving that 
G is not planar). By flipping pairs if necessary we assume Sj, is 
empty. Moreover, only one of Sj +1 and Sj +1 contains back-edges 
whose lower incidence is greater than b (for otherwise we have a con- 
tradiction proving that G is not planar). Up to a flipping of the pair 
of stack, we may assume this is Sj. Then we fuse Sj, . . . , S® into Sj 
andCS( e ) = ((5 1 ,5 1 1 ),(S 2 ,5 2 1 ),...,(^°,5 J 1 )). 

— similarly, if CS(e 2 ) = ((!?, 7?), (T°, I*), . . . , (T°, 7; 1 )), let q = 2 if 
T® or T{ has a bottom edge whose lower incidence is a and let q = 1 
otherwise. Then it should be checked that for q < q' < p, one of 
S®, and Sq, is empty (for otherwise, we have found a contradiction 
in the constraints proving that G is not planar). By flipping pairs if 
necessary we assume S®, is empty. It should be noticed that if q = 2, 
then one of Tj 3 and T\ is empty. By flipping the pair if necessary, we 
may assume this is T-f. Then T\ contains exactly the back-edges in 
L(e 4 ). 

Then we fuse T^ 1 , . . . , T p x into Sj and, if q = 2, we add to the one of 
Si and S\ which is not empty the edges from T°. 

As the number of pairs of stacks decreases at each "fuse" , as no new pair 
of stacks is created and as the initial number of pairs of stacks is 0(m), 
the global time spent in this step by the algorithm is 0(m). 

• Deletion of the back-edges whose lower incidence is u: by monotonicity, 
we only have to check the top of the stacks in the top pair of stacks. As 
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every back-edge will be deleted exactly once by the algorithm, the global 
time spent for deletions is 0(m). 

As the stacks used by the algorithm have their elements linked in a top- 
bottom manner, it is easy to keep, while "deleting" back-edges, the link between 
edges that have to get the same A- value and to add a special type of links between 
some edges which have to get different A- value. This way, a spanning forest of 
the constraint graph is maintained with a constant time cost per deletion, which 
allows to propagate A-values after the planarity testing phase. A A-valuc for 
each back-edge being computed, the embedding follows and is easily computed 
in linear time. 
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4.3 Example 



We consider a sample planar graph with an arbitrary Tremaux tree (see Fig.[SJ). 
The right ol the figure display the information computed by the two first steps of 
the algorithm: the status of the edges (block/thin/thick) and a TT-precedence 
order -<*, represented here as a circular order of the outgoing edges. 




Figure 5: A sample graph with a Tremaux tree T (on the left). The 
block/thin/thick partition and a TT-precedence order -<* (on the right). 

The algorithm first computes the CS's of the back-edges. Then it proceeds 
by iterating on the vertices whose outgoing edges are all processed. At each 
such vertex v, the CS of the incoming tree edge is computed, using a Merge and 
possibly a Deletion step (denoted by — ►). In the following table, the CS's are 
represented as sequences of double stacks. 
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v = 9 



CS((8,9)) = 



v = 10 


Cb((8, 10)) = 


(10,6) 






(10,7) 




7J = 5 


CS((4,5)) = 


(5,1) | 




j 


(5,2) | 




= 8 


CS((7,8)) = 


(10,6) 
(9,6) 






(9,7) | 


(10,7) 
















(10,6) 
(9,6) 




(Deletion) 


v = 7 


CS((6,7)) = 


(7,3) 






(10,6) 
(9,6) 
















(7,3) 




(Deletion) 


« = 6 


CS((4,6)) = 


(6,3) 
(7,3) 








v = 4 


CS((3,4)) = 


(5,1) | 






(5,2) | 








(5,1) | 






(5,2) | 




u = 3 


CS((2,3)) = 


(5,1) | 






(5,2) | 








(5,1) 




(Deletion) 


v = 2 


CS((1,2)) = 


(5,1) 









(6,3) 




(7,3) 





(Deletion) 



(Deletion) 
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